BILINEAR R-PARITY VIOLATION IN RARE MESON DECAYS 
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Abstract. We discuss rare meson decays K'^ — > 7r~£+£'+ and D+ K~l'^l''^ 
(i, I' = e, /i) in a supersymmetric extension of the standard model with explicit 
breaking of _R-parity by bilinear Yukawa couplings in the superpotential. Estimates 
of the branching ratios for these decays are given. We also compare our numerical 
results with analogous ones previously obtained for two other mechanisms of lepton 
number violation: exchange by massive Majorana neutrinos and trilinear R-parity 
violation. 

In the standard model (SM), the lepton L and baryon B numbers are con- 
served to all orders of perturbation theory due to the accidental C/(1)l x C/(1)_b 
symmetry existing at the level of renormalizable operators. But the L and 
B nonconservation is a generic feature of various extensions of the SM [1]. 
That is why lepton-number (LN) violating processes have long been recog- 
nized as a sensitive tool to put theories beyond the SM to the test. One of 
the most well known process of such type is neutrinoless double beta decay 
(A, Z) — > (A, Z + 2) + + e~ that has been searched for many years (see, 
e.g., [2] and references therein). 

In Refs. [3,4] rare decays of the pseudoscalar mesons K, D, Dg, and B of 
the type 

M+ M'-£+£'+ {£, t' = e, /i) (1) 

mediated by light (m^y ^ Tni^mii ) and heavy (toat 3> ttim) Majorana neu- 
trinos were investigated. The indirect upper bounds on the branching ratios 
for the decays (1) have been derived taking into account the limits on lepton 
mixing and neutrino masses obtained from the precision electroweak measure- 
ments, neutrino oscillations, cosmological data, and searches of the neutrinoless 
double beta decay. These bounds are greatly more stringent than the direct 
experimental ones [5]. 

In Refs. [6, 7] we considered another mechanism of the AL = 2 decays 
(1) based on the minimal supersymmetric extension of the SM with explicit 
i?-parity violation (/iMSSM, for a review, see [8]). i?-parity is defined as 
R = (— where B, L, and S are the baryon and lepton numbers 
and the spin, respectively. The SM fields, including additional Higgs boson 
fields appearing in the extended models, have R = 1 while i? = — 1 for their 
superpartners. 
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The most general form for the -R-parity and lepton number violating part of 
the superpotential is given by [8] 

Wft = £a/3 (^IXijkLfL^E, + X'ij.LfQ'^Dk + e^L^H^^ . (2) 

Here i,j,k = 1.2,3 are generation indices, L and Q are SU{2) doublets of 
left-handed lepton and quark superfields {a,/3 = 1,2 are isospinor indices), 
E and D are singlets of right-handed superfields of leptons and down quarks, 
respectively; is a doublet Higgs superfield (with hypercharge Y = 1): \ijk{= 
—Xjik), X'ijk ^'^d ei are trilinear and bilinear Yukawa couplings, respectively. 

We assumed in [6, 7] that the bilinear couplings are absent at tree level 
(ei = in Eq. (2)). As well known they are generated by quantum corrections 
[8] but the dominant contribution of the tree-level trilinear couplings to the 
phenomenology is expected. In the present report, we investigate the case of 
tree-level bilinear couplings: 7^ with A = 0, A' = 0. For this case, trilinear 
couplings cannot be generated via radiative corrections. 

The bilinear terms in the superpotential (2) induce mixing between the SM 
leptons and the MSSM charginos and neutralinos xj'^ in the mass-cigenstate 
basis and lead to the following AL = ±1 lepton-quark operators [9, 10]: 

+f3^DkPRU^UL + liiiVkPRl^hi + P^uPne'^dL) + H.c. (3) 

Here Pl,r = (1 T7''')/2, "/^ = the constants k„ {n = 1,2,3,4), 

f^k' Pki (^'* = 1)2,3) and (i'^ depend on the elements of the mixing matrices 
diagonalizing the neutrino-neutralino and the charged lepton-chargino mass 
matrices. 

The Lagrangian describing the bilinear mechanism of the decays (1) is 

C = CLH + CsM+C~g+C<,. (4) 

In addition to the LN violating part (3), it includes the SM charged-current 
interactions 

CsM = -^W+^'(^9e{x)^^PLl{x) + Y,'i{x)li.PLVgq,q'{x)) + li.c., (5) 

t 9,9' 

where i = e, ij.,t; q = u, c, t; q' = d, s, b; Vqq' is the CKM matrix, and the MSSM 
gluino-quark-squark and neutralino-quark(lepton)-squark(slepton) interactions 

[11] 

Cg = -V2gs^-^iq„L~g^^m - q.Rg^^H'k) + H-c, (6) 



4 

£x = V2g^{eLn{ll')^LXni'L + CRnWi'RXni'R) + Hx. (7) 
n=l 

Here (Ar)° 6 are the 3x3 Gell-Mann matrices (r = 1, . . . , 8) with color indices 
a,b= 1, 2, 3; the neutraUno coupUng constants are defined as 

ELnW = -T3(^)7V„2 + tan^iy(T3(V') - Q(V'))iVni, Efln W = tan^tyAT^i, 

where Q{4') and T3{tp) are the electric charge and the third component of the 
weak isospin for the quark (lepton) field ip, respectively, and Nnm is the 4x4 
neutralino mixing matrix. 

The leading order amplitude of the decay (1) is described by 9 diagrams 
shown in Fig. 1. The hadronic parts of the decay amplitude are calculated 




Figure 1; Feynman diagrams for the decay M+ M'~ + + in the bilinear |?MSSM. 
Bold vertices correspond to Bethe-Salpeter amplitudes for mesons. There are also crossed 
diagrams with interchanged lepton lines. 

with the use of a simple model for the Bethe-Salpetcr (BS) amplitudes for 
mesons as bound states of a quark and an antiquark [12] (sec also [3,4]). In 
addition, taking into account that the meson mass niM <C mw TmsusY , where 
mw is the W boson mass and msuSY ^ 100 GeV is the common mass scale of 
superpartners, we neglect momentum dependence in the propagators of heavy 
particles (sec Fig. 1) and use the cfFcctivc low-energy current-current inter- 
action. In this approximation the decay amplitude does not depend on the 
specific form of the BS amplitude and is expressed through the known decay 
constants of the initial and final mesons, /m and 



Finally, for the total width of the decay (1) we have obtained 



Tw = r(M+ ^ M'-l+l'+) = (1 - -dee' V 



n=l 
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Here 



= y ^^^2 [l_(/i^ + /i_)(2^)-l]2 ^l_(/^ + /_)(2^)-l] 

1+ 

x[ih+-z)ih.-z)il+-z)il.-z)f' (9) 

is the reduced phase space integral with h± = {1 ± mM' /itim)^ and l± = 
[{me ± to^/)/tom]^; A^c = 3 is the number of colors. 

For the numerical estimates of the branching ratios (BRs) , B^/ = r«' /Ftotai i 
we have used the known values for the SM couplings g and gs, meson decay 
constants, meson and Icpton masses [5], and a typical set of supersymmctric 
parameters [13]: a) the MSSM parameters: TOq = 70 GeV, /i = 500 GeV, 



M2 = 200 GeV, tan (3 = 4; b) the RPV parameters: |A| = {Yfi^i IMl^V^^ 



£2 



£3; 



0.1 GeV^ with lOAi = A2 = A3, |e| = ELi = |A| with ei 
M3 = {g1 / g'^)M2 at the clcctrowcak scale. Using the MSSM mass formulas [11] 
with the gluino mass rUg = M3, the masses of squarks and neutralinos and the 
elements of the neutralino mixing matrix were calculated numerically for the 
above set of parameters. The results of the calculations with the use of Eq. (8) 
are shown in the fourth column of Table 1. In the second and third columns 
of the table, the present direct experimental bounds on the BRs [5] and the 
indirect bounds for the Majorana neutrino mechanism of the rare decays [4] 
are shown, respectively. We see that the BRs for the bilinear RPV mechanism 
are much smaller than these upper bounds. For comparison, the trilinear RPV 
mechanism leads to the upper limits on the BRs of order 10^23 (10-24) for K 
[D) rare decays with the use of conservative bounds on the trilinear couplings 
\K3kK'3'k' \ ~ 10"^ But for more stringent bounds |A'A'| < SxlQ-^ [8,14], 
B«'(tri^MSSM) < IQ-^s. 



Table 1: The branching ratios B , for the rare meson decays M+ — > M' 



Rare decay 


Exp. upper 
bound on B«/ 


Ind. bound 
on B«/ {umSM) 


(bi^JMSSM) 




6.4 X lO-i'^ 


5.9 X 10"^^ 


3.6 X 10-'*'' 


7r^/i+//+ 


3.0 X 10-" 


1.1 X 10-^'' 


1.0 X 10-4" 




.").0 X 10"-^" 


5.1 X 10--^ 


4.2 X 10-^" 


D+ K-e+e+ 


4.5 X lO""^ 


1.5 X 10"^^ 


1.6 X lO-''^ 


D+ K-fi+n+ 


1.3 X 10"'' 


8.9 X 10-^4 


1.5 X lO"** 


D+ K-e+n+ 


1.3 X 10-4 


2.1 X 10"^=* 


3.1 X 10-4** 
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